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Abstract 

X 

We perform a general study of primordial scalar non-Gaussianities in single field inflationary 
models in Einstein gravity. We consider models where the inflaton Lagrangian is an arbitrary 
function of the scalar field and its first derivative, and the sound speed is arbitrary. We 
find that under reasonable assumptions, the non-Gaussianity is completely determined by 5 
parameters. In special limits of the parameter space, one finds distinctive "shapes" of the 
non-Gaussianity. In models with a small sound speed, several of these shapes would become 
potentially observable in the near future. Different limits of our formulae recover various 
previously known results. 
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1 Introduction 



Inflation has been a very successful paradigm for understanding otherwise puzzling aspects 
of big bang cosmology [1,2]. It can naturally solve the flatness, homogeneity and monopole 
problems that otherwise seem to require a very high degree of fine tuning for the initial state 
of our universe. Furthermore, inflation generically predicts almost scale invariant Gaussian 
density perturbations [3] , consistent with experimental observations of the Cosmic Microwave 
Background. Future experiments can constrain and distinguish between inflationary models 
in several ways. 

In the next few years, we can expect to see increasingly precise determination of the 
scalar spectral index n s and its running^]. Planck will also lower bounds on the tensor-to- 
scalar ratio r well into the regime favored by models of large- field inflation [6]. However, 
n s and r are just two numbers. Although their precise determination will be a tremendous 
achievement, it will leave considerable ambiguity in reconstructing the correct inflationary 
model. 

In contrast, the non-Gaussian component of the scalar fluctuations is characterized by 
a three-point function which is, a priori, a nontrivial function of three variables (momen- 
tum magnitude and ratios) on the sky. Furthermore, as demonstrated conclusively in [7,8], 
slow-roll models where the density perturbations are produced by fluctuations of the infla- 
ton itself, predict negligible non-Gaussianity. A detection of non-Gaussianity by the next 
generation of experiments would therefore strongly favor either an exotic inflationary model, 
or a model where density perturbations are generated by other dynamics (as in curvaton [9] 
and modulated reheating [10] scenarios). A crude measure of non-Gaussianity is the number 
/jvl- Values of \/nl\ > 5 would almost certainly indicate some novelty in the dynamics of 
the inflaton itself. 

In this paper, we determine the most general non-Gaussian perturbations possible in 
single-field inflationary theories. We assume that the inflaton itself generates the density 
perturbations, and that the Lagrangian is a function of the inflaton and its first derivative 
alone. Under these assumptions, we prove that the full non-Gaussianity (at the first order 
in various slow-variation parameters) is actually captured by five numbers. These numbers 
characterize the three-point function of fluctuations of the inflaton (or more precisely, its 
gauge-invariant analogue). For models with c s « 1, which are known to produce the 
most significant non-Gaussianities, the result is stronger: the leading non-Gaussianity is 
characterized by two numbers, and two different possible qualitative shapes in momentum 
space; the subleading non-Gaussianity is characterized by three more numbers. 

The elegant, gauge invariant calculation of the non-Gaussianities for slow-roll models ap- 
peared first in [7,8]. The result is that slow-roll models produce a primordial Jnl of O(10~ 2 ), 
too small to measure. Creminelli stressed that models where the effects of higher derivatives 
are important may give larger f^L (under the assumption that the non-Gaussianity is not 
diluted by the inflationary expansion itself, as happens in some models [11]). However, this 

1 Already, there are strong hints that models with a red spectrum are preferred [4]. Moreover, if the 
spectral index runs from blue to red, then there should be an approximate coincidence [5] between the 
length scale k at which n(k) — 1 = and the length scale at which the tensor to scalar ratio reaches a 
minimum. Such a coincidence of scales, if observed, can put constraints on inflationary model building. 



2 



effect can only reach /tvl ~ in the regime where effective field theory applies [12]. A 
proof of principle that significantly larger f NL can occur in sensible models was provided by 
the work of Alishahiha, Silverstein and Tong [13], who found that a fairly concrete string con- 
struction [14] could yield substantially larger f NL » 1. Different models in this class were 
constructed [15, 16] in the context of warped compactification. The large non-Gaussianities 
in these models are compatible with the current observational bound, but potentially ob- 
servable in future experiments [13,17]. (Another interesting model also providing large Jnl, 
but so far resisting embedding in a UV complete theory, appears in [18]. This model does 
not belong to the general class that we will study in Einstein gravity.) Our results build on 
these papers and further significant work by Seery and Lidsey [19], who found the fluctua- 
tion Lagrangian to cubic order for a general class of Lagrangians, and of Babich, Creminelli 
and Zaldarriaga [20], who emphasized the importance of analyzing the full shape of non- 
Gaussianities in /c-space. 

There are several motivations for completing such a general analysis. Firstly, it provides 
a null hypothesis against which to compare any future measurement of the non-Gaussianity. 
Secondly, several string-inspired models, if realized in nature, will give rise to a very char- 
acteristic measurable non-Gaussianity. String theory is relevant here because models with 
significant non-Gaussianity tend to be governed by higher-derivative terms, and a UV com- 
pletion is needed to make sense of such models (indeed, this is true even of slow-roll models 
with negligible non-Gaussianities [21], since the slow-roll conditions are sensitive to Planck- 
suppressed corrections to inflaton dynamics). Our analysis should make it straightforward 
to work out predictions for any such models. When we wish to provide specific examples, 
we use DBI inflation [14,15] and K-inflation [22]. 

In addition, non-Gaussian fluctuations could contain a signature of any departure of the 
inflaton from its standard Bunch-Davies vacuum. This has been suggested as a possible 
signature of trans-Planckian physics, and there has been much debate of the plausibility of 
such modifications; some representative references are [11,23-32]. It is a simple matter to 
translate any modification of the inflaton wavefunction into a modification of the three-point 
function, so the non-Gaussianities could serve as a test of any proposed modification. 

Finally, the structure of the three-point function can in principle be determined by 
dS/CFT [33] or its generalization appropriate to models with c s << 1. This could pro- 
vide a useful laboratory for studying holographic descriptions of dS space. This perspective 
could be useful even if there is no exact relation between the dS gravity theory and a dual 
field theory, since the useful aspects of the duality for this purpose are purely kinematical. 
Recent work in this direction appears in [34,35]. Our work on this connection will eventually 
appear in a companion paper [36]. 

The organization of this paper is as follows. In §2, we present the general class of 
Lagrangians that we will analyze (those that are an arbitrary function of a single scalar 
field and its first derivative), and define the notation we will use in the rest of the paper. 
While our analysis applies much more broadly, in §3 we describe three well-studied classes of 
inflationary theories whose non-Gaussianities we shall discuss in detail as special examples: 
slow-roll models, DBI models, and power-law K-inflation models. In §4, we find the cubic 
fluctuation Lagrangian in appropriate gauge-invariant variables for the most general single- 
field Lagrangian, and compute the non-Gaussianities. Our main result is that there are only 
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a few basic shapes, governed by 5 parameters in the most general model. We evaluate our 
results for the three special examples in §5, and present the different qualitative shapes of 
the non-Gaussianities that may occur. In §6, the effects of putting the inflaton in a vacuum 
other than the Bunch-Davies vacuum are described. We conclude in §7. The reader who 
is interested only in the class of Lagrangians studied and the general structure of the non- 
Gaussianity for this class, can confine her attention to sections §2 and §4 (which are more 
or less self-contained). 



2 Inflation models with a general Lagrangian 

To set up our notation, let us first review the formalism in [37] where a general Lagrangian 
for the inflaton field is considered. The Lagrangian is of the general form 

S = I J d'x^-g [M 2 pl R + 2P(X, 0)] , (2.1) 

where is the inflaton field and X = —^g^d^d^cf). The reduced Planck mass is M pi = 
(8nG)~^ and the signature of the metric is (—1, 1, 1, 1). The energy of the inflaton field is 

E = 2XP X - P , (2.2) 

where P t x denote the derivative with respect to X. Suppose the universe is homogeneous 
with a Friedmann-Robertson- Walker metric 

ds 2 = -dt 2 + a 2 {t)dx 2 3 . (2.3) 

Here a(t) is the scale factor and H = ^ is the Hubble parameter of the universe. The 
equations of motion of the gravitational dynamics are the Friedmann equation and the 
continuity equation 

SM^H 2 = E , (2.4) 
E = -3H(E + P) . (2.5) 

It is useful to define the "speed of sound" c s as 

2 dP P x 

Cs= dE = P~x + 2XP XX (2,6) 
and some "slow variation parameters" as in standard slow roll inflation 

H XP X 



e = 



H 2 M 2 t H 2 ' 

e 



c s H 



(2.7) 
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These parameters are more general than the usual slow roll parameters (which are defined 
through properties of a flat potential, assuming canonical kinetic terms), and in general de- 
pend on derivative terms as well as the potential. For example, in DBI inflation the potential 
can be steep, and kinetically driven inflation can occur even in absence of a potential. We 
also note that the smallness of the parameters e, rj, s does not imply that the rolling of 
inflaton is slow. 

The primordial power spectrum is derived for this general Lagrangian in [37] 

P C = ^ E^ = ^H 1 

k 36tt 2 M*c s (P + E) %n 2 M 2 pl c s e ' 1 ' ; 

where the expression is evaluated at the time of horizon exit at c s k = aH. The spectral 
index is 

dhxPjj , , 

n -- 1 = liu± = -*-•>-»• 

In order to have an almost scale invariant power spectrum, we need to require the 3 pa- 
rameters e, rj, s to be very small, which we will denote simply as 0(e). We note that in 
inflationary models with standard kinetic terms the speed of sound is c s — 1, but here we 
do not require c s to be close to 1. For example, in the case of DBI inflation, the speed of 
sound can be very small. In the case of arbitrary c s , the formula (12.81) (12. 9i) for the power 
spectrum and its index at leading order is still valid as long as the variation of the sound 
speed is slow, namely s<l. We will discuss this in more detail in Sec. 14.11 

The tensor perturbation spectrum Pj? and the tensor spectral index ht are given by 

ph 

and they satisfy a generalized consistency relation = — 8c s riT- This is phenomenologically 
different from standard inflation when the speed of sound is not one. 

3 Several classes of models 

In this section, we review three types of single field inflationary models. We discuss the basic 
setups and results of the corresponding effective field theories. These models will be used as 
primary examples after we work out the general expression for non-Gaussianities. 

3.1 Slow-roll inflation 

Slow-roll inflation models are the most popular models studied in the literature. The effective 
action takes the canonical non-relativistic form 

P(X,0) = X-^(0) . (3.1) 
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One achieves inflation by starting the inflaton on top of a flat potential V ((/)). The flatness 
of this potential is characterized by the slow roll parameters 

Ml (V- 2 

Vv = M p 2 ^, (3.2) 
which are required to be much less than one. The energy 

E = X + V^V (3.3) 

is dominated by the potential and the sound speed c s = 1. During inflation the inflaton 
speed is determined by the attractor solution 

V 

(3.4) 



3H 



This condition relates the slow roll parameters in (13.21) to the slow variation parameters in 

e = e v , r} — -2rj V + Ae v . (3.5) 

The primordial scalar and gravitational wave power spectrum are both determined by the 
potential 

p c = 1 y3 (ze) 

k " 12tt 2 M6 V' 2 ' 1 ' 
h 2V 

Pk = 3^iWj- (3 ' 7) 
The spectral indices and the running can be computed using the relation 

dink = Hdt = ^d<\> , (3.8) 



and we get 



dink ~ pl \ V 1 
din P f h 

dink " ~ 1v1 p 1 V^ ' 



V' 2 


V"\ 




+ 2 — 




v ) 


V' 2 V" 


vv 




V 2 


V' 2 





"r = -n-f = -Ml— . (3.9) 



There has been significant effort invested in developing slow-roll models of inflation in 
string theory. Some fairly recent reviews with further references are [38,39]. 
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3.2 DBI inflation 



DBI inflation [13-16] is motivated by brane inflationary models [21,40-43] in warped com- 
pactifications [44-49]. In particular, strongly warped regions or "warped throats" with ex- 
ponential warp factors, can arise when there are fluxes supported on cycles localized in small 
regions of the compactified space. A prototypical example of such a strongly warped throat 
is the warped deformed conifold [50,51]. The effective field theory of compact models con- 
taining such throats [48,49] has been explored in detail in [52]. In the slow-roll paradigm, in- 
flation can happen when a brane is approaching anti-branes in a throat if the potential is flat 
enough. However, this is non-generic [21]. Both the degree of tuning involved, and various 
possible ways of engineering flat potentials, have been discussed in the literature [21,53-59]. 

Perhaps the most interesting idea, which relies upon dynamics distinct from the usual 
slow-roll paradigm, arises in the DBI model. In this model, the warped space slows down the 
rolling of the inflaton on even a steep potential. (This "slowing down" can also be understood 
as arising due to interactions between the inflaton and the strongly coupled large-N dual 
field theory). This scenario can naturally arise in warped string compactifications [15]. The 
inflaton </> is the position of a D-brane moving in a warped throat. In the region where the 
back-reaction [14, 16, 60] and stringy physics [16, 17] can be ignored, the effective action has 
the following form 

S=^j d'x^R- J d 4 x v ^[/(0)- 1 V / l + /(0)^0^0-/(0)- 1 + ^(0)] . (3.10) 

The above expression applies for D3-branes in a warped background where f(<p) is the 
warping factor. We will first express the results in terms of a general f(4>). For an AdS-like 
throat, f(<f) — (where A in specific string constructions is a parameter which depends on 
the flux numbers) Two situations have been considered in the literature: 

• In the UV model [13, 14], the inflaton moves from the UV side of the warped space to 
the IR side under the potential 

V{<f>) ~ ^m 2 2 , m>M p //VX. (3.11) 

In this case the inflaton starts far away from the origin and rolls relativistically to the 
minimum of potential at the origin. 

• In the IR model [15, 16], the inflaton moves from the IR side of the warped space to 
the UV side under the potential 

V{<f>) ~ V - ^m 2 <j) 2 , m~H. (3.12) 

The inflaton starts near the origin and rolls relativistically away from it. 

The evolution of the inflaton in both cases was studied and the resulting power spectra 
were computed in [13-16]. Stages of DBI and slow-roll inflation can also be smoothly 
connected to each other [15,62]. 

2 This is a good approximation if we assume that the last 60 e- foldings of inflation occur far from the tip 
of the throat. Otherwise, inflationary observables may depend on the details of the warp factor [61]. 
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In the following we summarize the basic results of DBI inflation, following [13-16], and 
using the general formalism developed in [37]. For the zero mode evolution, we assume the 
inflaton (p is spatially homogeneous and denote X = <p 2 /2. The pressure P and energy E are 



p = -/(^rvi - 2xW) + f^r 1 - v{<j>) , 

E = 2XP, X ~P=—r ===; - + V{4>) , (3.13) 

V 1 - 2a/(0) 

and the speed of sound c s 

c s = y/l-fifM) • (3-14) 

In DBI inflation, the scalar rolls relativistically and a speed limit can be inferred by 
requiring positivity of the argument of the square root in the DBI action. So in this limit 
c s < 1, we can approximate the inflaton speed during inflation by 

±-^= . (3.15) 



^7W) Vx ' 

It is easy to see that the requirement e C 1 (or equivalently \E + P\/E <C 1) implies that 
the potential energy V(<p) dominates throughout inflation despite the fact that <fi is rolling 
relativistically. Hence, the Friedmann equation and the continuity equation reduce to 

V'W = -3H^=, (3.16) 

where we have used the universal speed limit relation ( 13.15P in the continuity equation. The 
number of e-foldings is computed as the following 



The scalar power spectrum and gravitational power spectrum are computed in the general 
formalism to be 

,c i e 2 _ /(0M0) 2 



k 



36tt 2 M* c s (P + E) 36vr 2 M p 4 z 



h 2E 2V(<j>) 

Svr 2 ]^ 37r 2 M* ■ 

The spectral indices and the running can also be computed using (I3.8p . 

cttntj _ V3M pl( j) 2 ( 4 | 2V 



(3.18) 



rflnfc y/\V \ (f) V 
dn s _ 3M 2 2 f 4 , 8(f)V' + 20 2 V" 30 2 V" 2 



rflnfc A V V V 2 V 2 V 3 



d\nP k h /3M 2 0V 

" T = 7inT = V^r^r' (3 - 19) 
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where we have evaluated f(<p) = A/0 4 . Using the equations of motion it is easy to verify 
that the gravitational wave spectral index satisfies the generalized consistency constraint 
$ = -8c s n T in [37]. 

We make two remarks regarding the result (13.191) . Firstly, as pointed out in [13], for 
the UV model, both the variation in the speed of sound and the small parameters e, rj 
contribute to the scalar spectral index and their effects cancel each other in the case of a 
quadratic potential V(4>) = ^m 2 (j) 2 , so that the spectral index is a second order quantity 
0(e 2 ) in this case. Indeed, we can directly see the cancellation from the first formula in 
( I3.19P for a quadratic potential. Secondly, in the IR model the potential remains V(<fi) ~ V 
during inflation, so the above expressions for the number of e-foldings, the scalar spectral 
index and its running ( 13. 19}) become simplified 



n s -l 



*f d4 1 [XVn I IXVn 



M pl \ 3 Mprt 



VXV v </>/ N t 

-2 
pi 



e 



dn s 3MU? 4, 4 



dink A v Vq N 2 



(3.20) 



In this IR model the gravitational wave production is very much suppressed compared to 

ph 

the UV model. This suppression is due to the consistency relation = — 8c s nr, and to the 

fact that the gravitational wave spectral index rij- is much smaller than the scalar spectral 
index n s — 1 in this case since |^r| <C |^|. 

A concern in DBI inflation is how to get the large background charge A (~ 10 14 ) which 
is needed to fit the field theory result to the observed density perturbations. Since this re- 
quirement just arises from requiring the compactification scale in the throat to be ~ Mqut 
(combined with the standard AdS/CFT relation between g s N and the compactification vol- 
ume), it seems very likely that model building could significantly reduce the apparent tune. 
For discussions of this issue, see [13,16,17]. 



3.3 Kinetically driven inflation 

One simple class of models which can give rise to large non-Gaussianities is the models of 
K-inflation, where the dynamics of inflation is governed by (non-standard) inflaton kinetic 
terms [22,37]. The Lagrangians giving rise to K-inflation are not radiatively stable, and so 
this mechanism is UV sensitive. There are as yet no convincing limits of string theory which 
give rise to K-inflation, but because the models are so simple, we analyze them in detail 
nonetheless. It would be very interesting to find controlled limits of string theory which give 
rise to such models. 

The simplest class of K-inflation models are the models of "power-law K-inflation." The 
Lagrangian for power-law K-inflation is of the form 

P(X, 0) = \^^^{-X + X 2 ) , (3.21) 
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where 7 is a constant, not to be confused with the Lorentz factor (l/c s ) that often appears 
in the literature of DBI inflation. (The form of the Lagrangian and our discussion can be 
straightforwardly generalized to a more arbitrary form where P oc f(X)/<p 2 [22].) Before 
describing the physics which follows from (I3.2ip . we should discuss some general concerns 
about K-inflation. The most obvious (also mentioned above) is that a Lagrangian of the 
form (13.211) is not radiatively stable, since it is not protected by any symmetry. (A shift 
symmetry of <fi could protect a Lagrangian of the form P(X) with generic coefficients). A 
second concern is that a reasonable exit mechanism for inflation must be provided. A third, 
related concern is that the dominant energy condition 

OP dP 

^>0, X--P>0 (3.22) 

is not satisfied by (13 . 2 1 [) for small values of X. Therefore, while in the inflating solution we 
will see that ( 13.221) is satisfied, one must provide an exit mechanism that changes the form 
of P drastically enough that the physics around flat space is sensible. We shall discuss these 
issues further after summarizing the key properties of the solution of interest. 
One solution to the equations of motion [22] is to take 

X = X = (3.23) 

which gives rise to an FRW cosmology with 

a(t) ~ t£ (3.24) 



for any < 7 < 2/3. The speed of sound following from (13.211) is 

cj = ^. (3.25) 

We are most interested in the regime with c s << 1. Therefore, we will focus on models with 
small 7, and sometimes expand formulae around 7 — > 0. 

3.3.1 The effective theory governing small fluctuations 

To get some intuition for these models, it is useful to construct an effective theory describing 
small fluctuations around the background inflating solution. The equations (I3.23f) imply that 

<p(t) ~ (1 + 7/8)* (3.26) 

where we have absorbed an overall constant into the definition of t, and have only written 
the solution to O '(7). Let us cast the Lagrangian into a more familiar form by performing 
the field redefinition 

$ = log(0) (3.27) 
valid for t > 0. Then the mini-superspace Lagrangian takes the form 

C = f{l)(-\& + \&eA (3.28) 
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/(T) ~ 7^ • (3-29) 



where f(p/) is the complicated 7-dependent prefactor in (j3.2ip ; for small 7 

16 

97 1 

Defining Y = —^g^d^d^, the full Lagrangian is just 

c = f(i) (-y + yV*) . (3.30) 

Now, we introduce the effective field it which describes small fluctuations around the 
solution via 

&(x,t)=ir(x,t) + $ (t) (3.31) 
where $0 characterizes the inflationary solution, with 

$ (t) = log ((1 + (3.32) 



(and therefore Yq = We will also find it useful to define 

1 

2 



Z = —g^d^irduit . (3.33) 



Expanding £ around this solution, we find a Lagrangian for 7r of the form 

£(7r,7r,V7r) = /( 7 )(£o + A + A) • (3.34) 

Here, the subscript on the £ denotes its order in the fluctuation tt. The first order term 
is guaranteed to vanish in the background $0, since it solves the equations of motion. To 
compute the Lagrangian up to third order in ir, we need the expansions of Y and e 2 *. Using 
the explicit solution (13.321) we see that 



and 



e 2# = t 2 (1 + 2) h + 2tx + 2vr 2 + |tt 3 H J . (3.36) 

A, = - J) , (3.37) 

A = (1 + ^ 2 + 2 Z + ^(1 + I) + 2^(1 + I) , (3.38) 

£3 = (1 + ^) + 2— + 27rZ + 27T7T 2 + 2t7T^ . (3.39) 

A \6 t z t J 

Here, the field 7r has dimension one, and all terms should be rendered dimension four by 
appropriate powers of M p i. One can use the relation H ~ ^ to replace explicit powers of t 
above with powers of H and 7, in estimating the size of various terms. 



We then find 



and 
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3.3.2 Basic phenomenology 



Here, we describe the rough phenomenology of a 'realistic' model of power law K-inflation. 
The exit from inflaton will be discussed in a later subsection. 

The power spectrum of these models was derived by Garriga and Mukhanov [37]. Working 
in the limit of small 7, their answer (equation (42) of [37]) becomes 

, c _ 1 2 GHl f k V 37 _ 1 2 Hi ( k ^ ~ 37 



Pk c s 3 7 vr UJ c s 3 7 87r 2 M^ UJ (3 ' 40) 

where H\ is taken to be the Hubble scale at the time of horizon exit for the perturbations 
currently at our horizon, and k\ is the the associated comoving wavenumber. 
It follows from (13^01) that the tilt 

n s -l = -3 7 + --- (3.41) 

which allows us to fix 7 ~ 1/60 using the central value of the spectral index in the WMAP 
results [4]. This justifies our use of perturbation theory in 7 in earlier equations. 

Using (13.25)) and (13.401) . as well as the fact that data determines ~ 10~ 9 at horizon 
crossing, we find 

3 

at horizon crossing, where M p i is the reduced Planck mass. Plugging in 7 ~ ^ as determined 
from n s , we see that the Hubble scale when the 60th from the last e- folding leaves our horizon 
is roughly H ~ 10~ 5 M P ;. So primordial gravitational waves will be unobservable in this 
model. 

The reader may wonder about the following. In these models, with c s << 1, the sound 
horizon where fluctuations freeze out at c s /H can be much smaller than 1/H. What happens 
if % < l p , i.e. H > c s M p il This would seem to give rise to a "trans-Planckian problem." 

However, it is easy to see that this regime cannot be reached in any reliable fashion. The 
power spectrum (13.401) makes it clear that 



H 2 ~ _L_ 7 3/287r 2 A*J x 10- 9 (3.42) 



Sp 1 / H 



p 7 s/4 y M} 



v 



(3.43) 



Hence, for H > j^-^Mpi, — > 0(1) and there is no good semi-classical description of any 
resulting region which exits from inflation. 



3.3.3 Exiting from K-inflation 

2 

Using a(t) ~ £ 3 t , the requirement that one gets 60 e-foldings starting from some initial time 
ti is simply 

\og{t f /U) = 90 7 • (3.44) 

Then, we need to arrange for an appropriate exit mechanism to kick in at the t j we determine 
in this way. 
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Here, we briefly describe a simple mechanism to exit from K-inflation, modeled on hybrid 
inflation [63]. Going back to the Lagrangian (I3.30p . we would like to arrange so that after 60 
e-foldings, at some specific value of 4>, the inflationary stage ends and the exit to standard 
radiation domination occurs. An easy way to do this, while fixing the problem that (I3.30p 
violates the DEC around X = 0, is to consider a more elaborate theory including also a 
second scalar field if). Then a Lagrangian of the schematic form 

C = -Y + F V* + (<9^) 2 + ($ 2 - $> 2 + ±-i?Y + a^ 2 - pMlf + ■■■ (3.45) 

can do the trick. Here M* is some UV scale, perhaps the string scale or the Planck scale. 
For small a/3 2 , assuming ~ M*, then right around the time when ($) ~ $* the ip field 
becomes tachyonic and condenses. It rolls to a vev (-0) ~ /3M*, and for (3 > 0(1), can 
correct the sign of the <fi kinetic term. 

The Lagrangian (13.451) can then support an early phase of K-inflation, and exit to a 
phase with normal kinetic terms for the various fields. The DEC is satisfied during both the 
inflationary phase and around the flat-space vacuum with (ip) = /3M*. It is an interesting 
question to check whether it is satisfied all along the trajectory from the inflationary phase 
to the final vacuum. 

4 Non-Gaussian Perturbations 

Now in the general setup of Sec. [2J we consider the non-Gaussian perturbations in the 
primordial power spectrum. There is a large literature on this subject, see e.g. [7, 8, 13, 
17,19,64,69,70]. However, most of the literature has been focused on the case where the 
speed of sound c s is very close to one, where the primordial non-Gaussianities are generally 
too small to be detected in future experiments. In addition, in some of the literature, only 
the perturbations in the matter Lagrangian are considered, but not the gauge invariant 
perturbation that remains exactly constant after horizon exitfj Here we will consider a 
general Lagrangian of the form (12.11) . and we will allow the speed of sound c s to assume 
arbitrary values, only requiring the parameters in (12. 7p (and one more parameter to be 
defined later) to be small and of order 0(e). We will calculate the three-point correlation 
function for the gauge invariant scalar perturbation ( following the approach of Maldacena 
[7]- 

It is useful to define two parameters following [19] 

E = XP X + 2X 2 P tXX = ^ , (4.1) 

A = X 2 P, XX + ^X 3 P, XXX . (4.2) 

In the case of inflationary models with P } x<f> — 0, the parameter A can be written in terms of 
the speed of sound and small parameters e, 77, s in (12.71) . However, for inflation models where 

3 For models with significant non-Gaussianity, this may be a reasonable approximation, since the contri- 
butions from the gravitational sector yield an /jv l which is too small to measure in any case. 
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P,x<j> 7^ 0, such as DBI inflation or K-inflation, there is no simple formula for the parameter 
A in terms of the slow variation parameters, and we must treat each model individually. 

To compute the Einstein action to the third order, it is useful to work in the ADM metric 
formalism 

ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) . (4.3) 

This formalism is convenient because the equations of motion for the variables A?" and N l are 
quite easy to solve. We will work in a comoving gauge where the three dimensional metric 
hij takes the form 

hij = a 2 e 2 % , (4.4) 

where we have neglected the tensor perturbations, a is the scale factor of the universe and 
( is the scalar perturbation, and remains constant outside the horizon in this gauge. The 
index on N l can be lowered by the 3-dim metric h^. The inflaton fluctuation 5(f) vanishes 
in this gauge, which makes the computations simpler. Using the ADM metric ansatz the 
action becomes 

S=^J dtd 3 xVh~N(R {3) + 2P) + - J dtd 3 xVhN- 1 (E ij E ij - E 2 ) (4.5) 

where we have set the reduced Planck mass M p i = 1 for convenience. The three-dimensional 
Ricci curvature is computed from the metric h^. The symmetric tensor Eij is defined 

as 

Ei j = \(h lJ -V l N J -V 1 N l ) . (4.6) 

The equations of motion for A" and N l are 

#(3) + 2 P - AXP tX - N- 2 {E ij E ij — E 2 ) = , 

Vj{N- x El) - V l (N- 1 E) = . (4.7) 

We follow [7] and decompose A^ into two parts Ni = Ni + d(i\} where d^N 1 = 0, and 
expand A^ and N l in powers of ( 

N = 1 + at + ot 2 + ■ ■ ■ , 
Ni = A/f } +NP + -- - , 

if> = A + i>2 + ■ ■ ■ , (4-8) 

where a n , N^ n \ijj n ~ 0(( n ). One can plug the power expansion into the equations of motion 
(14.71) for A^ and N l . At first order in (, the solutions [7, 19] are 



«i = 4, A/f } =0, Vi = -4+X> d 2 X = a 2 - 2 C, (4.9) 



s 



after choosing proper boundary conditions. 

In order to compute the effective action to order 0(( 3 ), as pointed out in [7], in the ADM 
formalism one only needs to consider the perturbations of A^ and A^ 1 to the first order 0(C)- 
This is because their perturbations at order 0(( 3 ) such as a% will multiply the constraint 
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equation at the zeroth order O(( ) which vanishes, and the second order perturbations such 
as «2 will multiply a factor which vanishes by the first order solution (14. 9p . So the solution 
(14.91) is enough for our purpose. This general conclusion in the ADM formalism can be seen 
as follows. 

What we have done so far is solve the constraint equations for the Lagrange multipliers 
N and Ni, which result from the variation of the action with respect to them 



5S 



d 4 x 5C(diN, N) 



d x 



_d(d l N) ON 



(4.10) 



where for simplicity we schematically denote N as either N or Ni. We expand N = N^ + 
AN = JVW + + + • • ■ for 



dC 



dC 



d(diN) 
dC 
dN 



dC 



|o + 
d 2 C 



d 2 C 



9(W U ' d{diN)d{djN) \° d i AN + 



dN^ + on 2 



(4.11) 



where the subscript means AN = 0. In order to get N^ we can neglect the terms involving 
d(d n~)8n • ^ n * s * s because this term starts from 0(Q as we can see from (14.51) and (j4.6j) . so 
it does not contribute to AiV at the first order 0(C)- 



The C(C°) terms in (ODjl are 

dC 

d x 



d(diN) 



dC 

\ , c od i 5N+— \ ^o5N 



. 



(4.12) 



This equation is consistent with the background equation of motion (12 .4p . The 0(C) terms 
in (14. 10p are 



d x 



dC 



dim) 



lo,C 



+ 



d 2 C 

dC . d 2 C . AT(1) 



dN 



dN 2 



5N 



diSN 
= , 



(4.13) 



where the subscripts C or C° denote the order of the perturbation C that we take. After 
integration by parts this gives the constraint equations for AiV at order 0(C), namely N^ . 
A similar procedure can be used to solve for AiV to order 0(C n ), namely up to N^ n \ 

We will next substitute these solutions for the Lagrange multipliers into the action and 
expand to order 0(C n ), where n > 3. We show that to do this the knowledge up to 

N (n-2) 

is enough. Let us look at the terms that possibly contain N^ n ~^ and N^ n \ 



AS 



d x 



dC 



d(diN) 



\ d i AN + - 



d 2 c 



2 d(diN)d(djN) 

BC 1 8 2 C 

^loAiV + ^44r| (AiV) 2 



oidiAN^djAN) 



dN' 



2 (dN)< 



(4.14) 
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The terms involving d ^J^ 9N = 0(C) are n °t written because such terms will not contain 
N(n-i) or N (n) at order 0(C n ). 

The following are all terms containing JVW (n > 2) in (I4.14p . 

/ A {a(iv)l»^■ iV, " , + §^ A^, " , ^ < 4 ' 15 ' 

Comparing with (I4.12p we know that this term vanishes. This is because after integration 
by parts, will multiply a term which is just the zeroth order constraint equation coming 
from (14.121) after integration by parts. Next we look at all terms containing N^ n ~^ (n > 3) 
in (jUlD , 

+ M l °^ iv< "" I, + l »^ iv< "" ,iv(11 } • (416) 

This term also vanishes because, after integration by parts, TV*™ -1 ) will multiply a term which 
is the first order constraint equation coming from (14.131) after integration by parts. 

Therefore our task is simplified. In order to expand the action (14.51) to quadratic and 
cubic order in the primordial scalar perturbation (, we only need to plug in the solution 
for the first order perturbation in iV and N l and do the expansion. The results can also be 
extracted^ from [7, 19] 

S 2 = [ dtd 3 x [a 3 4c 2 - ae(<9C) 2 ] , (4.17) 



c 2 



//•3 o 3 , 
dtd 3 x[-eaC(dC) 2 - a 3 (S + 2X)j^ + — pCC 2 

+ t^(3C - ^{did^did^ - d\d 2 ^) - 2a- 1 d i ipd l (d 2 ?P} , (4.18) 

where ( is the derivative with respect to t. One can decompose the perturbations into 
momentum modes using 

Uk = [ d 3 x C(t,x)e^ k - X . (4.19) 



4.1 The quadratic part 

To solve the quadratic part of the action ( 14.171) we define 



av^2e 

v k = zu k , z EE . (4.20) 



Note that there is a typo in term in Eq (44) of [19], it should be -^ijfd ipid ipi instead 
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This brings the equation of motion for the perturbation ( to a simple form 



z" 



v'i + cjk 2 v k -—v k = } (4.21) 

where the prime denotes the derivative with respect to the conformal time defined by dt = 
adr, r = -(aH)- 1 ^ + 0(e)). The leading order of z"/z = 2a 2 H 2 {\ + 0(e)) is contributed 
by the scale factor a which has the strongest time dependence. If the sound speed varies 
slowly enough, the leading behavior of the equation (14.211) is given by a Bessel function. We 
write it in terms of the Fourier modes of (, u k , using (I4.20p . 

u k = u(t, k) = lH (1 + %kc s r)e-^ . (4.22) 

Here we have made the approximation that the sound speed changes slowly, so this solu- 
tion has oscillatory behavior but with a frequency that is slowly changing due to the time 
dependence of c s . This requires 

- krAc s < c s kAr . (4.23) 

That is, the phase change in (14.221) caused by the change of c s is much slower than that 
caused by the change of r. This condition can be brought to the form 

- < H , (4.24) 

which is just the condition for small slow variation parameter s in (12. 7p . 

From Eq. (I4.22p we can see that, before horizon exit c s k > aH, u k is oscillating and its 
amplitude is decreasing proportionally to r. For v k , this is the leading behavior in flat space 
and we have chosen the standard Bunch-Davies vacuum. (We will discuss non-Gaussianities 
for other choices of vacua in Section [6]). After the horizon exit c s k < aH, u k remains 
constant [3]. This is most easily seen from (I4.2ip where we can neglect the second term, 
and we see that (for the growing mode) v k oc z k so that u k = constant. This conclusion is 
not going to be changed by the higher order interactions [7], because, as we can see from 
the interaction terms (14.181) . they involve either spatial derivatives, which can be neglected 
at super-horizon scales, or powers of time derivatives starting from second order. So after 
horizon exit the leading value of u k is determined by (I4.22p at r pa with the rest of the 
variables evaluated at c s k = aH. We emphasize here that the validity of our analysis only 
requires the variation of sound speed to be slow; the sound speed can be arbitrary. For our 
later purposes, the first order corrections to the leading behavior (I4.22p is also important. 
We work this out in Appendix [Al 

Now, we follow the standard technique in quantum field theory and write the operator 
in terms of creation and annihilation modes 

C(r, k) = u(t, k)a(k) + u*(t, -k)a f (-k) (4.25) 

with the canonical commutation relation [a(k), a^(k')] = (2n) 3 5^ (k — k'). 
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4.2 The cubic part 

The cubic effective action in (I4.18P looks like order O(e ) in the slow variation parameters. 
In slow- roll inflation, as emphasized and demonstrated in Ref. [7], one can perform a lot of 
integrations by parts and cancel terms of order O(e ) and 0(e). The resulting cubic action is 
actually of leading order 0(e 2 ) in slow roll parameters. A similar analysis can be performed 
for the general Lagrangian in Ref. [19], as well as in the case of interest here where the sound 
speed is arbitrary. Except for terms that are proportional to 1 — c 2 or A, the rest of the 
terms can be cancelled to the second order 0(e 2 ), 

S 3 = [dtd 3 x{-a\X(l-±)+2\)^ + ^(e-3 + 3c 2 s )U 2 

+ °±{e-2s + l- c 2 s )((d0 2 - 2a- 2 C,{dQ{dx) 
ci ci 



where the variable \ is defined in Eq. (14.91) . and in the last term 



SL . fdd 2 X 



5( 1 V dt 



+ Hd 2 x-ed 2 () , (4.27) 



/(C) = ^C 2 + ^CC + ^kpi-( d O(dO +d- 2 (d i dM(d j O)] 

+ Awl(90(dx) ~ d-^d&id^x))] . (4.28) 

Here d~ 2 is the inverse Laplacian, f^|i is the variation of the quadratic action with respect 
to the perturbation £, therefore the last term which is proportional to can be absorbed 
by a field redefinition of (. It can be easily shown that the field redefinition that absorbs 
this term is 

C - Cn + /(CO ■ (4-29) 

For the correlation function only the first term in (I4.28P contributes since all other terms 
involve at least one derivative of £ that vanish outside the horizon. The three-point function 
after field redefinition ( — > ( n + -^z( 2 becomes 

(C(X 1 )C(X 2 )C(X 3 )) = (Cn(xi)Cn(x 2 )Cn(x 3 )) 



+ ^((C„(x 1 )Cn(x 2 ))(Cn(x 1 )Cn(x 3 )) + Sym) + 0( V 2 (P<) 3 ) .(4.30) 

We proceed to calculate with the above cubic terms (14.261) . The terms in the last line 
in ( I4.26P are all of subleading order in slow variation parameters (12.71) . The interaction 
Hamiltonian from the leading terms to 0(e 2 ) is 

H mt (t) = - [ rf 3 x{-a 3 (S(l - i) + 2A)^ + ^(e - 3 + 3c^)CC 2 

+ ^(e-28 + 1- cMd() 2 - 2a^((d0(dx)} ■ (4.31) 
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One then computes the vacuum expectation value of the three point function in the inter- 
action picture that characterizes the primordial non-Gaussianities 

(C(t,k 1 )C(t,k 2 )C(t,k 3 )) = -i / < ^([C(t,k 1 )C(t,k 2 )C(t,k 3 ),^(t / )]) , (4.32) 

J t 

where to is some very early time when the vacuum fluctuation of the inflaton is well within 
the horizon, and t is a time about several e- foldings after the horizon exit. Translated 
to the conformal time r = —l/(aH), we can in a good approximation take the integral 
over conformal time r from — oo to 0. We follow the standard procedure and compute the 
contributions from various terms. In the following we first evaluate the leading contributions 
of each term. 

1. Contribution from ( 3 term. We denote K = k\ + k 2 + fc 3 , and find 

2Xc 2 H 2 e f° adr 



i(c 2 s - 1 + ^)^«(0, k!)«(0, k 2 )«(0, k 3 ) / 



H 3 

x uu y, ^jau y, ^jau ^ (27r) 3j3 ( y kQ + C.C. 

(IT (IT (IT < J 



du*(r, ki) du*(r, k 2 ) du*(r, k 
dr dr dr 

gl(cj - 1 + ^(arWk, + k 2 + k 3 )<n ■ (4.33) 



i=l 

L 2 



2. Contribution from term. We find 

i4(e-3 + 3c^)w(0,k 1 )w(0,k 2 )w(0,k 3 ) / a 2 rfr 



r.y */ , ^du*(r, k 2 ) rfu*(r, k 3 ) Wri s 3t . 3/ v^, \ 

x 2(u*(r, k x ) ^-^ + sym)(27r) 3 5 3 (^k i ) + c.c. 



1G( v v£ " 3 + 3c^)(2 7 r) 3 5 3 (k 1 + k 2 + k 3 )(Q p) 

s i=l * 

x (^ + ^ + sym). (4.34) 



X K 2 
3. Contribution from C(<9C) 2 term 

H 4 



, 6r , c4 (e - 2« + 1 - c») (2vr) 3 5 3 (k 1 + k 2 + k 3 )(H ^) 



i=i 



//, , w v . hk 2 + hk 3 + k 2 k 3 . hk 2 k 3 
x ((k 1 -k 2 )(-X+ + __) + sym ). (4.35) 



4. Contribution from ((d()(dx) term 

H 



i=i 



k x ■ k 2 )A| fci + fc : 



x ( v " 1 ^ (2 + ^p) + sym). (4.36) 
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5. Contribution from field redefinition £ — > £ n + j^Cn- 



4 C 2 S n - 

pi 



I il^ (27r)353(ki + k2 + k3)i m + sym) • (4 - 37) 

As a first step, we add all these leading contributions together. After some simplification, 
we find 

(C(kx)C(k 2 )C(k 3 )) = (27r) 7 <5 3 (k 1 + k 2 + k3)(Pf) 2 ^3 A, (4.38) 
where the above contributions to A are organized as follows 

+ (H(4E^+^E^+gE^ 

x * ' \ i>j i^j i 

+ £ f 4 E - r< E + E ■ (« 9 ) 

4.3 Correction terms 

We notice that, in (14.391) for general c s , the first two lines are not of the same order of 
magnitude as the last three lines. The former are 0(1) while the latter are 0(e). So for 
(? s -C 1, it is clear that the first two terms dominate. If for completeness one is interested in 
the full result to 0(e), however, small corrections to the first two lines may compete with the 
last three lines. This means that one must calculate the sub leading terms (of order 0(e)) 
for the first three integrations in Sec. 14.21 

In obtaining (I4.39j) . we treat all the slow- varying parameters in the integrand as constant, 
and we use the leading order solution (I4.22p . The corrections come from several sources. 

Firstly, there are corrections to the leading order u(t, k) in (I4.22p . which we work out in 
Appendix [A], 

where y — ^fM. 
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Secondly, various parameter in this solution as well as others in the integrand, including 
H, c s , A, e, are all time-dependent subject to the slow variation conditions. We Taylor 
expand such functions as 

f(r) = f(t K ) + ^(t-t K ) + 0(e 2 f) 

= fir K )-%^-^- + 0{e'f) , (4.41) 
at ti K tk 

where the reference point tk is chosen to be the moment when the wave-number K = 
k\ + &2 + &3 exits the horizon. 

Thirdly, the scale factor a also receives 0(e) correction, 

1 ' +-^—Ht/t k ) + 0(€ 2 ) . (4.42) 



H K T H K T H K T 

We then consider all types of corrections in the first three integrations (14. 33ft . (14.341) and 
(14.351) in Sec. 14.21 We leave the details of these calculations to Appendix [B], and summarize 
the final results in the following subsection. 

4.4 Summary of final results 

To first order in the slow variation parameters 0(e), the three-point correlation function of 
the gauge invariant scalar perturbation ( for a general single field inflation model is given 
by the following: 

(C(k 1 )C(k 2 )C(k 3 )) = (27r) 7 5 3 (k 1 + k 2 + k 3 )(Pi) 2 ^3 

x (A x + A c + A + A e + A v + A s ) , (4.43) 
where we have decomposed the shape of the three point function into six parts 

a, = g-i-Ap-p-^MMM, (4 . 44 ) 



A = (i-l-^J (eFx € + r]Fxr, + sFxs) 



(4.45) 



A' 



(4.46) 



+ (eF ct + V F CV + sF cs ) , 

V i *>3 f 

A v = ^I> 3 ) ' ( 4 - 48 ) 

A, = sF s . (4.49) 



21 



The definitions of the sound speed c s , £ and A are 

2 _ P,x 



c 



s Px + 2XP XX ' 

2 i 



E = IPx + 2A^P XX , 

A = X 2 P XX + ^X 3 P, XXX . (4.50) 

The definitions of the four slow variation parameters are 

H e c s , A , A _ „ , 

e = — — , 77 = — , s = , 1 = t-=. (4.51) 



is defined as 



H 2 ' ' e# ' c,.ff ' APT 



Hi 



P^ = -^T^^- . (4.52) 

Note that H, c s , e, A and £ in this final result are evaluated at the moment tk = — K l +0(e) 
when the wave number K = ki + k 2 + k 3 exits the horizon Kc s k = clrHk-i as indicated by 
the subscript K. Various P's are functions of ki, whose detailed forms are given in Appendix 

EH 

These results clearly illustrate that very significant non-Gaussianities /jvl ^> 1 will arise 
most easily in models with c s < 1 or | > 1 during inflation (while conventional slow- 
roll models enjoy c s — 1 and A/S = 0). We also see that for this wide class of models, 
the functional form of the leading non-Gaussianity is completely determined in terms of 5 
numbers: c s , ^, and the three slow variation parameters e, r\ and s. Note that in (I4.39P we 
are assuming that these parameters do not vary over the few e-foldings we see close to the 
horizon. If they do, the parameter counting becomes a little more complicated, and one can 
in principle extract more information from these results (by studying the running of the non- 
Gaussianities). Indeed after Taylor-expanding some slow- varying functions in working out 
all the correction terms in Sec. 14.31 one more parameter I comes up. However as illustrated 
in (14.441) . it happens that one can absorb I in ^ as ^[2 — (3 — 2ci)Z], where Ci is the Euler 
constant. (Later we will often refer to it as ^ for simplicity.) The error introduced to (jl? 



after this absorption is of order 0(e 2 ). So to the first order 0(e) that we are interested, our 
final result is still parameterized by five numbers. 



5 Size, shape and running of the non-Gaussianities 

In the previous section we have obtained the most general form of the primordial three-point 
scalar non-Gaussianities up to first order in slow variation parameters f)4.5ip in single field 
inflationary models, where the matter Lagrangian is an arbitrary function of the inflaton 
and its first derivative. This non-Gaussianity is controlled by five parameters — three small 
parameters e, rj and s, the sound speed c s and another parameter A/S. There is a large 
literature studying non-Gaussian features in models belonging to this class. One interesting 
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feature of our result is that we can take different limits and smoothly connect various previous 
results. We will also explore regions which have not been studied before. Before giving several 
major examples, we first discuss some general features of the non-Gaussianity obtained in 
Sec. El 

The correlation function in Sec. 14.41 is a function of three momenta forming a triangle. 
Therefore generally there are three interesting properties — the magnitude of the function, 
its dependence on the shape of the triangle and its dependence on the size of the triangle. 
Namely, these quantities determine the size, shape and running of the non-Gaussianity. 

To discuss whether a non-Gaussianity is large enough to be observed we first need to quote 
the experimental sensitivities. The non-Gaussianity of the CMB in the WMAP observations 
is analyzed by assuming the following ansatz for the scalar perturbation 

c = a - \f NL cl , (5.i) 

where (l is the linear Gaussian part the perturbations, and Jnl is an estimator parameter- 
izing the size of the non-GaussianityJ§ This assumption leads to the following three-point 
correlation function 

(CCkOCCkaXO*)) = (2 7 r) 7 5 3 (k 1 + k 2 + k 3 )(-^/^(P fc C ) 2 )^| ■ (5.2) 

Notice that this shape is different from any of the shapes in Sec. 14.41 except for A v . But we 
can set up a similar estimator Jnl for each of those different shapes of non-Gaussianities to 
parameterize its magnitude. This matching is conventionally done for the equilateral triangle 
case k\ = k 2 = k 3 . We then have 

fm = -^(3-l-|)+(3-2c0^, 



NL 



NL 



5 




~81 




35 / 




108 \ 






0(6) 





eX 



JnT = 0(e). (5.3) 

The bound on the parameter f NL from data analysis depends on the shape of the non- 
Gaussianities which we will discuss shortly. The current bound is roughly < 300 for 
the first three shapes [65], and |/jvl| < 100 for the rest [4]. A non-Gaussianity is potentially 
observable in future experiments if |/jvx| > 5 [66-68]. 

The magnitudes of A t , A v and A s are unobservably small, of order 0(e). The sizes 
of A are determined by e/c 2 s) rj/c 2 s and s/c 2 s . So in order to make the magnitude of these 
functions larger, we need the denominator c 2 s to be less than at least one of the slow variation 
parameters. (A similar conclusion for A/S term.) It is very interesting to construct such 
models, and we will show an example in this section. The most significant non-Gaussianities 
come from A\ and A c when c s C 1 and/or A/S 3> 1. 



3 The sign convention of Jml here follows Ref. [7], and it is opposite to the WMAP's sign convention. 
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As in Ref. [20], to show the shape of A(ki, k 2 , k 3 ), we present the 3-d plot x 2 x x 3 1 A(1, x 2 , x 3 ) 
as a function of x 2 = k 2 /ki and x 3 = k 3 /k\. The shapes of — A\/kik 2 k 3 and A c /kik 2 k 3 are 
shown in Fig. [1] and [2j We can see that A\ an d A c have overall similar shapes, but with 
opposite sign. The shape of — A t /kik 2 k 3 is shown in Fig. [3j The shapes of A v , A s are 
similar to that of A e up to a sign, in the sense that they all roughly approach to a pole in 
the squeezed limit, e.g. when k\ = k 2 and k 3 — > 0. 

For a small sound speed, both the leading non-Gaussianity A\, A c and the subleading 
non-Gaussianity A Q are potentially observable. Their shapes in the squeezed limit are all 
similar, 

Ay,o ^ h _ ^ /g ^ 

Details are shown in Appendix IB. 2[ 

The size of the non-Gaussianity also depends on the scale that we measure. Analogous 
to the spectral index, we define [17] 

_ d\n\f NL \ 

n NG -l = . (5.5) 

For example, if the main contribution to /tvl comes from a small sound speed, then n^G~ 1 ~ 
—2s. So in this case a measurement of the running of the non-Gaussianity directly tells us 
one of the slow variation parameter s. 

5.1 Slow-roll inflation 

We now reduce Eq. (14. 43ft to the slow-roll case. In this case the deviation of the sound speed 
from one is very small. We denote u — 1 — i < 1. Assuming u = 0(e), we can neglect 

both A, fl35S} since s ps ^ < 0(e) and *4 since .A = 0{ue). As in Ref. [19], if we further 
assume P % x<i> = 0, the relations 

X= e -(^(l-u)s- u y *x = -~ (5.6) 
follow. In this limit the Eq. (jP4l - ([4~4"9~|) become 
A = - I -^—s + 



+ 



+ r] 



\ i>j i^j % 

(4E*? + gE^4E% 2 ) 

\ i i^j i>j / 
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Figure 1: The shape of —A\/kik 2 k 3 
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Figure 3: The shape of — Ae/k^k^ 



Taking into account of the convention difference Ash = 4A urs, we recover the result Eq.(83) 
of Seery and Lidsey. 

Further setting the sound speed to be 1 (therefore u = and s = 0), using the relation 
(13.51) and taking into account of the convention difference ^Maidacena = 8^4 OU rs, we recover 
the result Eq.(4.6) of Maldacena. 

5.2 DBI inflation 

In Ref. [13,17,69], the three-point correlation functions in various models with potentially 
larger non-Gaussianities have been calculated by considering the perturbations in the matter 
Lagrangian only In general such a procedure is not valid because it is only the gauge 
invariant combination ( that will remain constant after the horizon exit. However (as realized 
by these authors) for the large non-Gaussianity case, due to large non-linear self coupling 
of the inflaton, the corrections coming from the gravitational part are relatively small, and 
this simpler method can be successfully used to obtain the leading behavior of the non- 
Gaussianities. Indeed, in this limit our leading behavior in A\ and A c , which is of order 
0(c~ 2 , A/S), recovers the leading behavior of the result in Eq.(6) of Gruzinov [69]. However 
this procedure does not guarantee that the subleading term in Ref. [69], which is of 0(1), 
will be correct. In fact, as we discussed at the beginning of this section, for < c 2 s < e, r], s, 
the subleading terms actually come from A and a subleading term in A\, which are of order 
cj 2 0(e,r], s) and should be observable. 

Now we consider the example of DBI inflation discussed in Sec. 13. 21 Interestingly for this 
type of inflation model, because of the Lagrangian (13.131) . the parameter A defined in (14. 2 p 



26 



is 2 

A = X 2 P, XX + ^X S P, XXX = ^(1 - c 2 s ) . (5.8) 

So the leading order contribution in A\ vanishes. The leading behavior of A c reproduces the 
result of Alishahiha, Silverstein and Tong [13,17]. Using the estimator defined in (15.31) . we 
have 

f° « 0.32^^ « 0.32c; 2 . (5.9) 

To estimate the subleading order corrections, let us look at both the UV and IR model. In 
the UV model, the sound speed is given by 

Pl\mM vi 




It is easy to see that whether the subleading order is of order 0(1) or 0(e/c 2 ) depends on 
the value of <fi. In the IR model, the sound speed is related to the number of e-foldings N e 
before the end of inflation by 

cj 1 « /3N e /3 , (5.11) 

where (3 parameterizes the steepness of the potential j3 = m 2 /H 2 , which is generically of order 
one as we know from the usual eta-problem in slow-roll inflation. Since e, rj, s = 0(N~ l ), 
we see that this model falls into the region where c 2 <C 0(e,r],s), where both the leading 
non-Gaussianity A c and the subleading A Q are observable. 

We can also compute the running of non-Gaussianities as considered in [17] 



n NG - 1 w -2s = -2\ — — ^— + -)■ 5.12 

\l f(cj>)v(4>)\2V v <p) y J 

Just like the spectral index, the running non-Gaussianities also provide a good probe of 
the effective potential V(0) and warp factor /(</>), and could (in an optimistic scenario) 
distinguish between different possible background geometries where the brane motion is 
occurring. 



5.3 Kinetic inflation 

We saw above that in one large class of models with measurable non-Gaussianities, the lead- 
ing order of A\ vanishes. However, generally it does not vanish for some other models such 
as K-inflation, and could in principle be comparable to the second term A c . An experimen- 
tal detection of the shape of large non-Gaussianities could therefore in principle distinguish 
between DBI inflation and more general models where the first term has significant contri- 
butions. 

Here, we first outline a quick estimate of the non-Gaussianities in K-inflation, and then 
describe the exact result following §4. 
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5.3.1 Crude estimate of non-Gaussianities 

Given the fluctuation Lagrangians (13.381) and (13.391) . we can do a simple estimate of the 
non-Gaussianity (following the general strategy also employed in [13, 18]). While this is 
not strictly necessary in view of the detailed formulae in §4, it is perhaps illuminating to 
understand in simple terms why these models have large /atl- 
The basic point is the following. We saw in (13.401) that 

p '~^© 2 ' (513) 

Using the fact that 

C ~ f 7T (5.14) 

and evaluating this on the inflationary solution, we find 



lM pl 

It follows that 



C ~ ~4r- ■ (5-15) 



(CO ~ ^(vnr) . (5.16) 

T 

Furthermore, given the overall factor of f(j) ~ in C, which in particular multiplies the n 
kinetic terms, one has the relation 

(tttt) ~ 7 2 (<5tt) 2 (5.17) 
where 5tt is a typical fluctuation of the n field during inflation. Therefore, estimating via 

P f ~ ^(tttt) ~ (5tt) 2 (5.18) 



7 2 



and using ~ -f^-, we see that 



3/2 

ftr ~ ^ ■ (5.19) 

Using the form of the modes one can also see that 

6-k ~ Hdir . (5.20) 

Now, we are interested in estimating the non-Gaussianity, say through a crude estimate 
of /nl- Since naively 

^ ~ fNL^ , (5.21) 

we can plug (15.191) and (15.201) into the fluctuation Lagrangians to estimate the non-Gaussianity. 
The term 2t-kZ in C 3 contributes fluctuations of size j^-^i^jz) 3 ~ ^jji (H 5 /M p i). The 

largest terms in C 2 scale like We then find 

Cs 1 H 



C 2 7 7/4 M pZ 
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(5.22) 



Because 

this translates into the rough estimate 

Jnl ~ - ~ \ ■ (5.24) 

This in fact reproduces the more detailed results of §4 when applied to K-inflation, though 
it does not give the (potentially very important) information about the detailed shape of the 
momentum dependence. 

5.3.2 Shape of Non-Gaussianities in K-inflation 

An order of magnitude estimate for the general size of the non-Gaussian signatures arising 
in K-inflation, appears in (15.241) . Here, we refine this estimate using the formulae of §4. 
The Lagrangian of power law K-inflation, at leading order in the expansion in 7, is given 

by 

1fi 1 

P(X,0) = — -(-X + X 2 ). (5.25) 



We then see that 



and 



32 I 2 8 1 

A " ~ ' (5 ' 27) 

In each case, the estimate after ~ follows from the fact that X = ^ + 0(7) on the inflationary 
solution. 

We find that, to leading order in 7, 



7 V K 3 

and 



Ax ^ l l(^M) (5.28) 



( k S k ' k * + 2IO ^ k * k * + 8 53 k i 

\ i>j iy£j i 



So unlike the other higher derivative model we've carefully examined, DBI inflation, these 
models receive a leading-order contribution from the shape A\. Also, in distinction to a 
general DBI inflation model in Sec. 15.21 the constant speed of sound implies that the non- 
Gaussianity does not run in power law K-inflation. 

Evaluating on equilateral triangles and comparing to the "local" form of non-Gaussianities, 
this translates to an estimated f^L of 

170 1 , , 

/m - - gj-- ■ (5.30) 
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This gives f^L ~ 125 for realistic models, which is allowed by current experimental bounds 
but would be easily detectable in future experiments. Comparing to the rough estimate 
(I5.24p . which was f^L ~ ^, we see that the two results agree up to an 0(1) coefficient 
(whose moderate smallness prevents the model from being excluded by current data). 

6 Non-Gaussianities as a probe of the inflationary vacuum 

There has been some interest in the question of whether we can observe trans-Planckian 
physics in the Cosmic Microwave Background radiation [11,23-32]. In this context, the 
assumption that the Bunch-Davies vacuum is the unique initial state of inflation has recently 
been questioned. While several plausible alternatives for the initial state have been suggested, 
its precise form is highly dependent on how we model the behavior of quantum fields at 
Planckian energies. What is universal however is that any deviation from the Bunch-Davies 
vacuum during inflation will result in modulations of the power spectrum [24,25,29], thus 
offering the exciting possibility of probing the initial state of the universe from cosmological 
measurements. 

We hasten to stress that the microphysics that determines the choice of inflationary 
vacuum is by no means understood. Here we put aside the conceptual issues associated 
with the choice of vacuum and its consistency, and simply approach the problem of vacuum 
ambiguity from a phenomenological perspective. To be specific, we explore the possibility 
of using primordial non-Gaussianities to test any deviation from the standard Bunch-Davies 
vacuum. As it turns out, the effect of deviation from the Bunch-Davies vacuum on the 
shape of the primordial non-Gaussianities is quite simple to compute within our formalism. 
A general vacuum state for the fluctuation of the inflaton field during inflation can be written 
as follows 

u k = u(r, k) = -^£= (C + (1 + ikc s r)e~ lkc ° T + C_(l - ikc s r)e lkc ° T ) (6.1) 

In the standard Bunch-Davies vacuum we have C + = 1 and C_ = 0. Now we allow a 
small deviation from the Bunch-Davies vacuum by turning on a small finite number C_, 
and calculate the corrections to the shape of non-Gaussianities we found in the Bunch- 
Davies vacuum. For simplicity we only consider the corrections to the leading order non- 
Gaussianities A\, A c in the small sound speed c s C 1 limit (the corrections to the other 
shapes of non-Gaussianities due to a non-standard choice of vacuum can be worked out 
by a similar procedure). The computations of the three point functions of the primordial 
perturbations essentially go through as before. The first sub-leading correction to Bunch- 
Davies vacuum result is simply to replace one of the three u(t, k)'s with its C_ component. 
(A correction of 0{CS) in u(0, k) gives a term which has the same shape as in the Bunch- 
Davies vacuum case. We do not include them here.) This does not change the common 
factor 3 A 3 but will simply replace one of the fc, in the shapes A\, A c with — fc,. We denote 

the corresponding corrections A\, A c . We immediately find the corrections as 

i o\ Qppp 

Ax = Jfe(ff-)(^-i-^F^ 

c\ S 2 
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Figure 4: The shape of \A\\/kik 2 k 3 



+ 77 u I u V* + 7 u I l I l vJ . ( 6 - 2 ) 



v (A:i + A; 2 - fc 3 ) 3 (h-k 2 + k 3 ) 3 (-h + k 2 + k 3 ) 3 

s p=l i>j i^j i 

We can estimate the size of the non-Gaussianities A\ and A c according to the WMAP 
ansatz. This estimate is usually done in the equilateral triangle limit; we find 

1 2A 
Inl = -5Re(C.)(- - 1 - -) , 

C s n 
25 1 

f c NL = T ife(C_)(--l). (6.4) 

If the sound speed c s is sufficiently small, the effects of a slight deviation from the Bunch- 
Davies vacuum is potentially observable by future experiments. We plot the shapes of the 
non-Gaussianities A\ and A c in Fig. H]and Fig. [5j We see the shapes of these corrections are 
very distinctive and in fact dramatically different from that of the DBI inflation or slow roll 
inflation. In particular, these shapes are highly peaked at the "folded triangle" limit where 
fc 3 w ki + k 2 for arbitrary values of k\ and k 2 - This feature is not shared by other known 
sources of non-Gaussianities, and so measurements of the shape of non-Gaussianities could 
in principle be an excellent probe of the choice of inflationary vacuum. 

Note that, while the rising behavior of the non-Gaussianity in the folded triangle limit 
is the signal of the non-Bunch-Davies vacuum, the divergence at the limit e.g. k\ + k 2 — 
&3 = is artificial. This divergence is present because we have assumed that such a non- 
standard vacuum existed in the infinite past. Realistically there should be a cutoff at a large 
momentum M for k/a, where A; is a typical value of fci,2,3- This amounts to a cutoff for r at 
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Figure 5: The shape of \A c \/kik 2 k3 

t c = —M/Hk. Since the integrand is regulated at r = —1/Kc s due to its rapid oscillation, 
if r c < —1/Kcg, the cutoff M has no effects to our calculation. That is, for K 3> kH/Mc s , 
we will see the behaviors shown in Fig. H] & [5] near the folded triangle limit. But within 
K < kH/Mc s , the cutoff takes effect first, the divergence behavior will be replaced. The 
details depend on the nature of the cutoff, e.g. a naive sharp cutoff will introduce oscillatory 
behavior. 

7 Conclusion 

The forthcoming suite of cosmological experiments will nail down with ever greater precision 
the parameters of the inflationary model that yielded our homogeneous, isotropic universe. 
Some measurements, like the value of the spectral index and the nature of its running, 
are guaranteed to occur. Others, like a detection of primordial gravitational waves, are 
not necessarily expected to occur on theoretical grounds (since models with very small r 
seem more natural as quantum field theories), but would be tremendously exciting and 
instructive if they do. The discovery of significant non-Gaussian scalar fluctuations falls into 
this latter category. While the simplest models of inflation do not produce this phenomenon, 
its discovery would tell us something qualitatively important about the inflationary epoch, 
and experiments sensitive enough to measure |/jvx| > 5 will be launched in the next two 
years. For this reason, we feel it is worthwhile to parametrize the reasonable possibilities, 
and understand the qualitative physics of the models that produce them. 

In this paper, we have taken some steps in this direction for generic single-field models. 
There are several clear directions for further work: 

• It would be nice to derive the same formulae governing non-Gaussianities as arising di- 
rectly from symmetry principles. Perhaps these would be encapsulated most neatly in a 
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hypothetical dual, non-gravitational theory. For the models with c s << 1, this theory may 
have novel properties. 

• Higher derivative terms play a significant role in the dynamics of those single-field models 
which produce striking non-Gaussian signatures. One class of models where such terms are 
important, the DBI inflation [13-16], has a reasonable microscopic justification in string 
theory. It would be interesting to find other examples where one can microphysically justify 
the study of dynamics that is very sensitive to higher derivative terms. 

• We have focused here on single-field models. It is a logical possibility that our 60 e- 
foldings arose from a multi-field inflationary model. This could be motivated if, for instance, 
r is measured to be non-negligible. In slow-roll models, measurable r implies an inflaton 
that traversed a super-Planckian distance in field space [71], as in chaotic inflation [72]. At 
least in string theory, this is difficult to accommodate in single-field models [73], but could 
conceivably happen in a multi-field setting [74]. For this and other reasons, it would be 
worthwhile to develop a general framework for analyzing non-Gaussianities in multi-field 
models. Examples of multi- field models with significant non-Gaussianity appear in [75-77]. 
Formalisms to compute the non-Gaussianities in large classes of such models are developed 
in [78-81]. 
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A Corrections to Uk 

In this appendix, we calculate the 0(e) correction to the solution of Eq 
the method of [82,83]. 
We define 

c s k 
V = ^H 

and write the equation of motion of the quadratic action 

z" 

V k + C s fc V V k = 



• f)4.2ip . generalizing 

(A.5) 
(A.6) 
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in terms of y. Note that generally c s is a (slowly varying) function of time. Using 



we get 



Z ^ = 2a 2 H 2 {l-\e + \^\s) + 0{e 2 ) (A.7) 



(1 - 2e - 2s) y 2 ^- S y^+ y 2 v k - (2 - e + ^ - 3s)v fc = . (Ai 



The solution of this differential equation is given by 

v k = y^ 1+s) [dH® ((1 + 6 + s)y) + C 2 H® ((1 + e + «)„)] , (A.9) 

where 

* = § +e + K- (A ' 10) 

The Bunch-Davies (BD) vacuum corresponds to C 2 = 0. To determine the coefficient C\, 
we need to look at the large k behavior of the equation (1A.6[) . 

v£ + c 2 k 2 v h = . (A. 11) 

The behavior is more general than the usual case where the sound speed is constant, because 
here we allow c s to vary slowly. Using a similar approach and defining y = —c s kr = 
(1 + e + 0(e 2 ))y, we get the solution for (IA.11I) with positive energy (BD vacuum), 

v k - * e W+-»-f4 , (A.12) 
y/2c s k 

up to a constant phase. Here the coefficient is determined by the quantization condition 
(Wronskian condition), v k ^- — Ufc-j^ = —i, to first order 0(e). Notice that in ( 1A.12I) . the 
sound speed c s runs as a function of y. So we can expand it as 



- e^-^, (A.13) 
where the subscript on c s denotes the evaluation at y . Expanding flA.9j) in the same limit, 

v k - dJl 71 1 y s/2 e'W^-Ml , y » 1, (A.14) 
V 7T VI + e + s 

we find 

C 1 = -^-l=(l+ e -+ 8 -)^K (A.15) 
2 Vc s fcK 2 2 

Note a convention for the variables used here: the variables such as c sk with the subscript 
k are evaluated at 

yo = = i ; (A.16) 
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in some later formulae, the variables without the subscript such as c s mean that they are 
functions of y. 

From the definition 

v k = zu k , z = , (A.17) 

we get the expression for u k to order 0(e), 

My) = -^^ s W^ l + h^ S ^ y3/2H " ] ((l + e + s)y) ■ (A ' l8) 

As an application, we use the y — > limit of Eq. (IA.18I) to derive an expression for the 
density perturbation to order 0(e). To do this, we use the expansion of the Hankel function 
in the i/<l limit, 

Hl l] (y) - -»• * , ! , ^ (Tf , (A.19) 



sin vtx T(— v + l) \2 
and get 

where 

c 2 = ci - 2 + ln2 « -0.73 
and Ci = 0.577- • • is the Euler constant. Hence the density perturbation is 



k 3 

'u k {y = Q)\ 



2tt 2 
I H k 



(l - (c 2 + l)e - |t7 - (| + l) a ) + G(e 2 ) . (A.2f 



This generalizes the result (31) of Ref. [82] to the case with running sound speed. 
B Details on the correction terms 

In this Appendix, we provide details on the correction terms in Sec. 14.31 Let us look at the 
first integration (14.331) in Sec. H~! 



where 



/d 
dr a A(r) JJw(0,ki)— M*(r,ki) ■ (2vr) 3 5 3 (^ k,,) + c.c. , (B.22) 
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We have evaluated the leading contribution of this integral in Sec. 14.21 The corrections come 
from several different places. 

The first is from the time variation in /i(r), 

/i(r) = h{t K ) + ^{t-t K ) + 0{e 2 h) 

= f 1 (r K )-^^-\n- + 0(e 2 f 1 ) . (B.24) 

Ot li K T K 

We choose to eventually evaluate all the variables at the time tk, which is defined as the 
moment when the wave-number K = k\ + hi + k% exits the horizon Kc s k = clkHk-, at which 

tk = + 0(e) . (B.25) 

Kc sK 

All the subscripts K denote the evaluation at the horizon exit point defined in flB.25j) . The 
df\ jdt can be expressed in terms of the slow variation parameters 

^ = ( V e + e 2 )(cj 2 - c; 4 ) + es(-2c; 2 + 4cJ 4 ) + 2/Atf- 2 + 6e\H~ 2 , (B.26) 
where we have defined 

(B.27) 

We assume that the time variation of A is slow, I = 0(e). Pluggingthe corrections terms of 
(1B.24j) into (1B.22j) and evaluating the rest in leading orders, we geto 

D ^(l-^c0((-e-,)(i-l) + Ki-2) + (2Z + 6 e )A^MM. (B .28) 

The second comes from the correction to the scale factor a ~ — j*-. This can be obtained 
from the relation 

dt 1 , /1\ 
dr = - = --d - (B.29) 



and the expansion 



H H K 

Integrating ( IB. 291) we can get the following expansion 

1 e e 



H 



' 1 + t{t - t K ) + 0(e 2 ) . (B.30) 



H K T H K T H K T 



+ -—\n(r/r K ) + 0(e 2 ). (B.31) 



6 The integration dxln(— x) e lx = ici — -| has been used. Similar types of integrations will be used 
later. 
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Plugging this correction term into (1B.22[) . we obtain 



w i )e (i-i-a 5MMM. (B , 2) 



2' V<* E/a- 2A ' 3 

The third comes from the correction term to w(r, kj) which we obtained in Sec. \M We 
first look at the corrections to the factor u(0, kj) in (IB.22|) . The corrections to the final result 
come not only from the corrections in the bracket of ( 1A.20j) , but also from the running from 
ki to K 

' 1 - (e+ - + -)ln-M + 0(e 2 ) . (B.33) 



^/C si €i y/CsK^K V 2 2 K 

So these add a correction term to the three-point correlation function^ 



A.A D (-3(c 2 + l)e-— ^-3(y + l)s - (e+ - + -)ln— — 
£7*. 2K 3 



x (B.34) 



We next look at the corrections to the factor jpu*(r, kj) in (1B.22H . To do this we use 
( 1A.18P and expand the pre-factor around k = K, 



2V2 ^e K c sK k^ 1 \ 2 2 2 2 t k J 

x y ^H^((l + e + s)y) . (B.35) 

Denoting Au(r, fcj) as the corrections to the leading order, we have 



Au*(r,k t ) = -I^=^4^-^) e - 
2 ^c sK e K k.' 

[—i(e + s) + (e + s)x + isx 2 
(He + - + -) - (e + - + -)x - «x s) In — 



7 To make sure that the expansion such as (|B.33[) is perturbative, we need fcj ^> C(i^e _1 / e ). Note that 
this condition still allows one of the momenta to be much smaller than the others, e.g. (kz/ki) 2 -C e. 

In order for the expansion such as (1B.31|) to be perturbative, we need rxe^ 1 ^ 3> r ^ r^e 1 ^. So it 
appears that the integration over r can only be taken from tkg 1 ^ to rjfe -1 ^. We first look at the upper 
bound. Since the mode ki exits the horizon at r, ~ — l/fcjC S j. At the upper bound of r, all modes have 
exited the horizon and their amplitudes are frozen. So the error introduced by including the integration 
from tke^ 1 ^ to is of order 0( ^ si K e -1 ^) ~ C(e _1 / e ). We next look at the lower bound. For the range 
of ki that we are interested in, at r ~ r^e 1 / 6 all modes are well within the horizon. Their contributions 
are regulated away due to their rapid oscillation. Therefore, to order 0(e^ 1 ^ e ), we can effectively take the 
integration range for r from — oo to 0. 
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where x = — kiC S KT. The above corrections include those in the first line of (IB. 351) . in y 3 ^ 2 
where y = —kiC S KT~(l — e — s In + 0(e 2 ) is used, and in H^*((l + e + s)y) which includes 
corrections in the index v and corrections in the variable y. Differentiate (IB. 361) . we have 



1 H 



K 



2 ^c sK e K 



3/2 



r(«-N)t. 



c S Re 



x - e + - + -) + (e + - + tex + sx 2 

[ v 2 T K 2 2x 

+ (ie H — ?7 is — sx)xln — 

2 2 tk 



rfz/ A 2 2 ;J 
The first two lines in the square bracket in (1B.37|) contribute 



(B.37) 



AA D 



3 / 1 



4 \c 2 s 



- 1 



2A 



K 



((3-6c 1 )e+(--3c 1 >?+(3c 1 



17, , k^k^k^ 



- — s 



K 3 



+ (<+H)(^EW-;|£W))- 

The last term in (IB. 371) involves special functions and contributes 



where 



AA D -ie 



Ri(fa, k 2 , h) 



V s ( 1 



2 2 y V c 



k 2 k 2 



- 1 



2A 



k 2 , h) + sym , 



K 



A U*( \f-l - k 2 + h , -i^+tAx! 

axi n [Xi)(l— i — Xi)e k i 



xHi%(x) + xHW = 2[v - \)H™ X 



(i) 



x 



dJ u {x) 
dv 

dN u (x] 
dv 

dJ v (x) 
dv 



dN v (x) 
dv 



-TTX 



[sina;Ci(2x) — cosa;Si(2a;) 



7ix I {cossCi(2x) + sina;[Si(2a;) — tt]} 



[cosa;Ci(2a;) + sin xSi(2a;)] , 



1 

-TTX 

2 



-TTX ) {sinxCi(2x) — cosx[Si(2x) — tt]} 



(B.38) 



(B.39) 



(B.40) 



Xi = -kic sK r , 

h(x) = -2ie ix + ie~ ix (l + ix) [Cx(2x) + iSi(2x)] - m sin x + in cos x . (B.41) 
We have used the following relations, 



(B.42) 
(B.43) 

(B.44) 

(B.45) 

(B.46) 
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The same procedure can be repeated for the second integration (I4.34p in Sec. 



-2i / dr a 2 / 2 (r) (u(0, k 1 )n(0, k 2 )u(0, k 3 ) 

J — oo 

w , , du*(r, ko) du*(r, kO , . ^^.v-v, , /T , 

x ^(r,^) ^-^ L^Z + sym ) ( 2 7r) 3 5 3 (^k,)+c.c. , (B.47) 



where 



/ 2 = -(3-3c*-e). (B.48) 
From the variation of /2, we get the correction term 

aa d - 3 S)( I - 1) - 1) f(i - 2 Cl) ^E^ - (» - <=i)^E*?«J 



(B.49) 



From the correction term to the scale factor a, we get 



From the correction term to u(0, fcj), we get 

AA D - 3(c 2 + l)e + — ±77 + 3(^ + 1)5 + (e + ^ + -) In " 



x 



4 V 2 ' v 2 y v 2 2' if 3 

(H.(££w-s*g««)- (R51) 



The correction to u*(r, fcj) in ( IB. 361) contributes 



and 



+ (-(1 + 2c x )e + (i - Cl )77 - (| + Cl )a)l £ fc?^ 

+ (cie-^l-cOTZ + ia + cO.)-^^^ 3 ) (B.52) 



^-HH^HKf 161 ^ 1 • (R53) 
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where 



d = Re 



k 2+ k 3 , 



dx\h*{xi)e % fe i 



(B.54) 



The correction to jpu*(r,ki) in (1B.37|) contributes 



AA D 



(± - l) (((2 + 4 Cl )e + (-1 + 2 Cl )r7 + (3 - 6 Cl )*)l £ k 2 k) 



3 / 1 
4 



+ (-2cie + (1 - aft + (-4 + 6c x ) S )^ J2 k i k . 



.3 



+ (2e + ?? + s)k 1 k 2 k 3 

+ (l-2c lW i5E^ + ^E^ 



fix 



A' " 



and 



where 



(xa- = -Kc sK t) 



A - 4D 3^- 1 ),( £ + i + 5) (A ' 1 + Sym) ' 



Mi = — fciRe / ax i — (k 2 + ki + ik 2 k 3 — xje fc i — 

Jn xi * k\ ax i 



(B.55) 



(B.56) 



(B.57) 



/o x i 

Notice that in (IB.55I) . the last term is divergent. This divergence is cancelled by the diver- 
gence that appears in (IB.57I) using the limit h(x) — > (—2 + ci)i + i ln2x + 0(x 2 ) as x — > 0. 
So we can re-define 



i- x e -ix K 



Mi + sym = Mi - Re(Af + A$) / (1% h sym 

Jo ^A 

to absorb this divergence. 

The following are the corrections to the third integration (14.351) 



(B.58) 



x u 



-2i / dr a 2 / 3 (r) ( w(0, ki)u(0, k 2 )«(0, k 3 ) 
*(r, ki) W *(r, k 2 )u*(r, k 3 )(-k 2 • k 3 ) + sym) • (27r) 3 <5 3 (]T k 4 ) + c.c. (B.59) 



with 



(B.60) 
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From the variation of /3, we get 

.A 3 ( ( 2-£ )(1 -^ ) + £ 



X 



V E *. 3 - 5*i w» + 5 E **s + E - V 1 E 

» i^jr' i>j 

(B.61) 



From the corrections to a, we get 



X 



( =r E *? - \ k **< + \ E mj - Hr 1 E *?*? + £ E 



From the corrections to u(0, hi), we get 

AA D -I(3(c 2 + l) £ +^ + 3(| + l) S + ( £ +| + £)ln^ 



X 



Corrections from u*(r, fcj) give 

x 



((3 Cl e + ^(-1 + Cl )77 + ^(1 + Cl )s) k ! 
\ - j 

+ (~3e- -?7- 2 S )5^^ 2 

3 5 

+ (3e + -rj + (- - 2c 1 )s)/c 1 A; 2 A;3 

+ (6(1 + 2 Cl )e - 3(1 - 2 Cl )rj + (5 + Q Cl )s)^J2 k i k J 

i>j 

+ (-6c ie + 3(1 - Cl )r) - (1 + 6 Cl )s)-^ £ fc?fcj 

- (i+ci)4e*.-- c ^e^) 



and 



4Vc? J K \ 2 2 
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(B.62) 



(| " l) (^E*? + |E *W " ^ E • (B-63) 

v s \ i i>j i+j / 



(B.64) 



AAD-U\-l) (e+^ + ^)(N 1 + sym) , (B.65) 



where 



iv ! = — > fc Re / ctxi^e fc i (—1 — z — ; xi H — ; — xf)/i (xi) . (B.66) 

2 ^ Jo x{ h fci 

B.l Final results 

Collecting all the results in Sec. 14.21 and this appendix, we get the final result: 
(C(k 1 )C(k 2 )C(k 3 )) = (27r) 7 5 3 (k 1 + k 2 + k 3 )(P^ 2 — 



ru 3 

x (A + Ac + Ao + A + A) + A) (B.67) 
where we have decomposed the shape of the three point function into six parts 

Ax = — + (3-2 Cl )Z- -hk 1 , B.68) 



A c = - - 1 



L ) (4 E + 2^ E W + g E *?) • (B.69) 
A = (\- l -^=r) (^Ae + vFx v + sFxs) 

+ {eF ct + V Fc V + sF cs ) (B.70) 

A = ^(gE^J > (B.72) 

A = sF s . (B.73) 
The definitions of the sound speed c s , £ and A are 

2 _ P,x 



A 



'2 ; 



£ = XP X + 2X^P XX , 

A ee X 2 P iXX + ^X 3 P tXXX . (B.74) 

The definitions of the four slow variation parameters are 

e "-^' ^df' z "apt- (R75) 



P^ is defined as 



pC s JL . (B.76) 
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Note that H, c s , e, A and £ in this final result are evaluated at the moment tk = — K x c +0(e) 

when the wave number K = ki + k 2 + k 3 exits the horizon Kc s k = clkHr. So p£ in ( 1B.76|) 
is defined differently from (1A.21|) . The various functions F are given by the following: 

3 9 39 3 hk 2 k 3 \ k\k\k\ 
— Ci — —c 2 -In 



1 2 z 4 2 *~ K 3 J K 3 

+ 4^ E W - ^ E W + **> to) , (B.77) 

i>j 

Fx, - ^-(3 Cl -f)M|M, (B.78) 
^ ss ±^ + (^-6)*^, (B.79) 



4 



hk 2 k 3 



2>K^ 1 3 v 2 2 z 1 3 y 2 L/ K 3 

{-\ E^ 3 + ^E k n - ^ E W) + 1^. **> **) . (b- 8 °) 



+ ( t + x E ^ + ( -f - j C2 + \>h E ^M-M^^P 
f-s E*. 3 + ^E - ^ E + w ■ < B - 81 ) 



+ In 



^3 



Fes EE "(— + -C 2 + -)Vfef--y^ + (£1 + 1)^^3 

V 16 16 16 ; ^ 1 16 ^ J V 4 16 ; 

+ { \ Cl + l C2 + z>ic E W + (f c i - i C2 " i } i E 

8 y J K ^ 1 8 K 2 ^ * J 



+ m 



f-^ E *? + ^ E - 4^2 E + **) , (B.82) 



^3 

F * = -7 c iE^ 3 + lE^^-I A; ^ 2A; 3 + (-2 + 2c 1 )-^^^ 



i^j i>j 
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+ (f-cO^M+(?-3co^f 



(B.83) 



where Ci = 0.577 • • • is the Euler constant and c 2 = c x — 2 + In 2 = —0.73 • • ■. The functions 
R(kx, k 2 , fa) and Q(fa, fa, fa) involve special functions, 

+ sym , (B.84) 

ih2+h3. xi dh*{xi] 



R(kx,fa,fa) 
Q(h, k 2 , fa) 



h 2 h 2 

— fciRe 

.3 , 7,3 > 



Re 



ax i \ 1 — i — x\ ] e 



,■ fc 2+ fc 3 



/o x l V "4 



- (k 3 2 + kl)Re 

h 2 h 2 

- -^Re 

+ sym , 



[°°dxi -ihi+H Xl ( .fa + fa k 2 k 3 2 \ 
/ — T e [-1 — 1 — - — Xi + —j-2-x]_\n 

Jo x i \ ki fa J 



x Xj 

(B.85) 



h(x) = -2ie ix + ie~ ix (l + ix)[Ci(2x) + iSi(2x)] - zvr sin x + ivrx cos a; , (B.86) 
Xi = -kiC sK T , x K = -Kc sK t . 

In all formulae, the "sym" stands for two other terms with cyclic permutation of the indices 
1, 2 and 3. 



B.2 The squeezed limit 

It is interesting to look at the behaviors of various functions in the squeezed limit (for 
example, k\ = k 2 and fa — > 0), because from them one can roughly know whether the shape 
of a non-Gaussianity is closer to the DBI type (Fig. |5J), or the slow-roll type (Fig. [3]). 

In slow-roll inflation, Maldacena has argued that the three-point function in the squeezed 
limit goes to [7] 



(C(k 1 )C(k 2 )C(k 3 )) - -(2vr) 7 ^3K - 1)^4 



(B.87) 



to all orders in the slow-roll parameters. This condition was generalized to general single-field 
inflation in Ref. [84]. 

Here we check this condition in the general single-field inflation model with large non- 
Gaussianities. For a large A/E 3> 1, the r.h.s. of the consistency condition starts from order 
(9(e) (does not include the terms P^P^), while the l.h.s. starts from 0(A/£) in A\ and 
(9(eA/S) in A a . So in order for the condition to hold, both the leading and subleading order 
terms in the three-point function have to vanish in the squeezed limit. It is not difficult to see 
that A\ vanishes in this limit. Interestingly, the subleading terms ( IB. 771) . ( IB. 781) and ( IB. 791) 
also vanish in this limit due to a cancellation from the special function R (see Appendix 

MM. 
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The case with a small c s <C 1 is similar. The condition (IB. 87ft requires the leading and 
subleading orders of the l.h.s. vanish. It is easy to see the leading order contribution A c 
satisfies this condition. For the subleading order 0(e/c 2 s ), (1B.80|) . (IB. 811) and (1B.82|) vanish 
in this limit due to a cancellation from the special function Q (see Appendix IB.3j) . Overall 

A Q goes as 



Ao h 
oc — . 

k\k 2 k 3 h 

A t and Ay are the same as the slow-roll case; in addition, since 



F s _ _^ 
k x k 2 k 3 Ak 3 



(B.88) 



^B.89) 



the order 0(e) terms on both sides of the condition also match, as we can see from fl2.9[) . 
So we have checked that all orders (the leading, subleading and next-to-subleading orders) 
of our full results satisfy the consistency condition. 

B.3 Some details 

In this section, we demonstrate some details that the first six F's vanish in the squeezed 
limit. To calculate the squeezed limit of the i?-term, we analytically continue the integrand 
in the convergent direction by x — > —ix and note the asymptotic behavior 



x 



Ci(-2ia:) - zSi(-2zx) tt e~ 2x (l + £>(-)) 

2 2x x 

So the e x terms cancel in the asymptotic behavior of h*(x), and we find 



OG 



h*(—ix) ~ xe x , x — > +oo . 
One can also compute the asymptotic of h*(x) near x ~ 0, and we find 

h*(x) = (-c 2 « -ilnx) + 0(x 2 ) . 



(B.90) 



(B.91) 



(B.92) 



We can now take a squeezed limit k 3 — > and k\ = k 2 = k in the i?-term. There are 3 terms 
in the symmetric rotation of the indices. Two terms are proportional to k^ and the integral 
is convergent because of the good asymptotic behavior (IB.91j) . (IB.92j) . so they vanish in the 
squeezed limit. The remaining term is 



R(ki,k 2 ,k 3 



—Re 

h 



dx-i 



2k 
h 



e'^ xa h*(x 3 ) 



[B.93) 



We only keep leading term in the squeezed limit and drop higher powers of k 3 . We can then 
change integration variable x 3 — > ||y and expand around k 3 = using the formula (IB. 921) . 
We find 



R(k 1 ,k 2 ,k 3 



Rc 



2k 



y + 0{kl) 



POO 

-i / dy (1 - y) e 
Jo 



y (c-zln(^)) 



+ 0(k 
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So the squeezed limit of the term F Xe vanishes 

,3 3 3 



g - I + ^)k 3 + 0{kl) = O(kl) . 



(B.95) 



The Q-term is more complicated. In the same limit 

Q 



+ 



2k 3 Re 
2k 3 



^f-2 2/ e^#^(|) + (2z-3,)e-- 



Re 



o y 

00 dxx 



o 



x 



dy 

Cl (l + ix!)/i*(xi) 



- -k 3 
3 



(B.96) 



where the 1st line comes from the first and second lines of (1B.85|) and their two cyclic 
permutations, the 2nd line comes from the last line of (1B.85j) and it cyclic permutation 
k\ — > k 2 , the 3rd line comes from the rest. Similar to the previous case, we have used the 
change of variable Xi and the expansion of the function h*(^j-y) around zero. To 

further evaluate (IB.96j) . we need the full expressions 



h*(x) 

dh* 
dx 



ie ix (l - ix)Ei(-2ix) 



2ie~ ix - ^(1 + ix)e~ ix 

(l--x- -)e- ix - ixe tx Ei(-2ix) . 
2 x 



Re i I ^r(l + ix)e- tx 



x- 



Plugging in these expressions, and making use of the integrals 

1 , 

dx Ei(— 2ix) = - , 
; 2 

2x) cos2x + (2 — nx) sin2x — si(2x) 



—i 



o 

dx 

~~2 



X 
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2 



we get 



So in the squeezed limit, 



Q -» —k* 
V 3 



IF. 



CI) 



2F r 



-5k 3 + -Q 



(B.97) 
(B.98) 

(B.99) 
(B.100) 
(B.101) 

(B.102) 
(B.103) 



Although not explicitly demonstrated, one expects that they vanish as O(k^) as we see 
explicitly in (IB. 951) . because in this limit the 0{k^) term vanishes due to isotropy [84]. 
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